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The paper deals with the relative convergence of a doubly-infinite set of equations per- 
taining to a boundary value problem in a waveguide. It is shown that the solution for the 
equations converges differently for different choices of the combination of the number of 
equations from the two sets. This is demonstrated by studying the convergence properties 
of the truncated set of equations as the number of equations is made indefinitely large. It 
is proven that the solution for the reflection coefficient becomes identical with the exact 
solution only with a particular choice of the ratio of the equations. This choice of the 
unique ratio is also shown to be consistent with the edge condition. 

1. Introduction 

In formulating electromagnetic boundary value problems we are often led to an infinite 
set of equations. In most cases, it is not possible to invert the infinite matrix and we are forced 
to resort to truncating the above to a finite size. We then solve the equations for a number of 
increasing sizes of the matrix and study the convergence of the solution. If at least the leading 
members of the unknown coefficients (which are usually of primary interest) tend to converge, 
we feel satisfied and assume that we have obtained a reasonably good approximation for the 
leading coefficients. 

It is the purpose of this paper to show in the first instance that for a particular doubly 
infinite set of equations associated with the bifurcation problem in a waveguide, there is a 
relative convergence of the solution, meaning that the solution converges to a different set of 
answers for every different choice of the ratio C=P/Q. P and Q are the numbers of equations 
from the first and second set, respectively, out of the doubly infinite set of equations. These 
equations for the bifurcation problem have been obtained by Hurd and Gruenberg [1954], 
who have also presented an exact solution of the infinite set through the use of calculus of 
residues. The above set, because it has a known solution, is particularly suited for out* purpose 
which is to demonstrate by considering several finite size matrices that there results a relative 
convergence when P-^ 00 and Q-^ 00 . Furthermore, we are also able to find the ratio of PjQ 
which yields the correct answer in the limit. 

In the second part of this paper we present a basis for choosing the correct ratio of PjQ 
when working with a truncated set. With this choice, when the size of the set is increased 
indefinitely while keeping the ratio C constant, the solution does converge to the correct answer. 
It is shown through the study of the asymptotic behavior of the higher order unknown coeffi- 
cients and the application of the edge condition, why only an unique choice of the ratio would 
make the solution asymptotically tend to the correct one and why otherwise an incorrect 
solution will result. 

2. Development of the Infinite Set of Equations for H-Plane Bifurcation 

The infinite set of equations which will be discussed here in connection with the problem 
of //-plane bifurdation in a rectangular waveguide has been derived by Hurd [1954]. We shall 
therefore skip the details and merely outline the procedure for their derivation. 
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Figuee 1. H-plane bifur -cation of a rectangular waveguide. 



The geometry under consideration is shown in figure 1. Assume that the incident wave 
from the negative ^-direction is a TE 10 mode with the electric vector parallel to the edge of the 
septum. It can be easily shown that the only nonzero field components are E y , H x , and H z and 
that they can all be expressed in terms of a scalar function <p=E y and its partial derivatives. 

The problem can be stated in terms of the equation 



V|rf>+*V=0 



(1) 



and the following boundary conditions on cp: 

(a) cp and W are finite everywhere in the region concerned, except at the edge of the bifur- 
cation at 2=0, x = b where V<p becomes infinite. 

(b) tp and V<p are continuous in the subregions and at 2=0. 

(c) <p behaves as an outgoing wave at large z apart from the incident field. 

(d) (p vanishes on the walls x=0, a and on the septum at x=b, z^>0. 

(e) <p satisfies the edge condition and hence goes to zero at the edge as r 1/2 , where r is the 
distance from the edge. V<p therefore goes to infinity as r~ 1/2 as r^O. 

It is easily verified that <p A , <Pb and <Pc expressed in the following equations satisfy (1) 
and the conditions (a), (c), and (d). The expressions are: 



(p A =Ae~ a i z sin (tx/cl)+^>2 A n e a n z sin (nrx/a) 



(Pb=^2 B n e &n z sin (nirx/b) 



(Pc=^i O n e y n z sin nir(x— b)/c 

71 = 1 



(2a) 
(2b) 
(2c) 



where 

A= amplitude of the incident field 

a n =[(nT/a) 2 -k 2 ) 1 ' 2 
P n =[(mr/by-JcT 2 
y n =[(nTlcy-k 2 Y' 2 
k = 2ir/\, X=free space wavelength. 

a, 0, and 7's are the mode propagation constants in the three regions. 

By applying the continuity conditions at 2=0 and subsequently equating the Fourier 
components of the resulting equations in the range 0-CxK.b and 6<Ca;<Cc, one arrives at the 
following doubly infinite set of equations after some manipulations (for details see Hurd [1954]): 
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s 



A n sin (nirb/a) A sin (71-6/a) 



a n —$ v &i+P P 



~ -4 n sin (rnrb/a) _A sin (7r6/&) 



2=1,. 



(3a) 

(3b) 



where 5^=0 for nj*l and equals 1 when n=l. 

The set of infinite equations (3) has been solved exactly by Hurd. We shall, however, 
concentrate on the solution of the truncated set of equations for various combinations of 
P and Q where these are the number of equations from the set (3a) and (3b) respectively. 
We proceed to do this in the following section. 



3. Solution of Truncated Set of Equations 

Consider the solution of the truncated set of equations 

^±Q A n sin (rnrb/a) A sin (wb/a) p 

h " <*n-P P «l + ft, P~h *,..•* 

r ^ A n sin (rnrb/a) _ A sin (wb/a) __ n 

£i _ «,- 7 , ^+7, 2-1,2, ••■V 



(4a) 
(4b) 



In attempting to solve the set of equations (4) we recognize first of all that the determinant 
of the equation is of a particular kind which is called a double alternant [Muir, I960]. Written 
explicitly, the determinant is 

1 1 1 1 



A= 



«i — ft" 



X P+Q 



-A 



a l — $P a 2 — @P 



a P+Q ftp 



«i--Ti 



«2— ?i" 



<*P+Q ^1 



«i— 'Vq 



X P+Q 



-y Q 



(5) 



Note that we have grouped the factors sin (rnrb/a) with the unknowns A n . 

It is observed that each element of the determinant is a reciprocal of the difference of two 
quantities, a r and f3 s or y t in general, only one of which, viz, the a, changes as one goes along the 
columns, whereas only /3 or y changes as one goes down the rows. Hence, if a r ->a n , where r 
and n are two different subscripts, then it is obvious that the determinant A— ^0 and that the 
zero is a simple one. This is effective to saying that (a r — a n ) for various combinations of r 
and n are the factors in the continued product expansion of the determinant. In a similar 
manner, we observe that (/3 r — /3 n ), (y r — In) and (f3 t —Vs) are also to be included in the product 
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expansion of the determinant. We have to make sure, however, that all such factors are in- 
cluded in the product and that none is repeated. It is also obvious that the denominator of this 
expression for A must contain the factors (a t —/3 s ) and (a t — y r ) for various combinations of the 
subscripts. The expression for the A developed on the basis of above arguments, is 

M, M-l P, P-l Q, Q-l P, Q 

n ia m -a n ) n (p m -Pn) n (y m -y n ) n (P m -y n ) 

. / -.\„ 771=71+1, 71 = 1 m = n+l, 72 = 1 771=71+1, 71 = 1 771 = 1, 71 = 1 / n\ 

&—V*-) M,P M, Q W 

n K-A) n (a m -y n ) 

771 = 1,71=1 771 = 1,71=1 

where M=P-\-Q; v is an even or odd integer, depending on P and Q, and as yet undetermined, 
although we shall not really need to find v. This is because in order to calculate the coefficients 
A T we need find only the ratio of the determinants A r /A where A T is the determinant obtained by 
replacing the rth column of A by the column representing the right-hand side of (4). It is 
fairly straightforward that 

A T sin (r7r&/a) =— — A(a r ^>—a 1 ) sin (wb/a) (7) 

where A(a r ->—a 1 ) is the determinant A with a T replaced by —a u We can therefore obtain 
through the use of (6) and (7) an expression for A r in the product form, which is 

A sin (nrb/a) _ . A(a r -^—^i) 
r sin (wb/a) A 

=-A n fo-*) n (7n ~ ar) M n K+C * l} (s) 

,1=1 (jffn+ax) 71=1 (7n + «l) 71 = 1 k-«r) 

where the superscript (r) on the last product implies that the factor corresponding to n=r is to 
be omitted. Equation (8) gives the desired solution for the truncated matrix. In particular, 
the expression for A 1 /A, which can be identified as the reflection coefficient R, is 

R= 4± = - n (ft *~ ai) n {yn ~ ai) n ((Xn+(Xl) 

A 72=1 (<Ui+Pn) n = \ («i+7n) 71=2 (« n — «l) 

____ £ (1-aSn) « (l-ajjn) ™ (l+«iAO (9) 

71 = 1 (l+«i/A) 71 = 1 (l+ai/7n) rc=2 (1—<*iA*») 

We also quote below for comparison and reference, the result arrived at by Hurd as a 
solution of the infinite set of equations. His expression has a sign error and corrected it reads 

B =- ^am-C:m-^l) X 6XP ^H-«+*ln (a/b)+cln (a/c)}] (10) 



where for instance 



n(co, a)= n (a p — oo)(a/p7r)e (ao)/Pir) 
p=i 



and the superscript (1) implies as before that the factor corresponding to p=l is to be omitted. 

4. Comparison of Solutions of Finite and Infinite Set of Equations 

In this section, we shall study the asymptotic behavior of the expression for the reflection 
coefficient R. We shall see that even when (P+Q) is very large, the expression (9) for the 
reflection coefficient yields different values for different ratios of P/Q. We shall then compare 
the answer for a particular choice of P/Q and see that it indeed converges to the exact value of 
R calculated through the use of (10). 
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Consider the particular case when all the mode propagation constants excepting a x are 
real and a x is purely imaginary. This is done merely for convenience, and the conclusions 
regarding the convergence phenomenon reached in this case will still be applicable to the 
general case when the above condition is not true. 

For this case, it is obvious from (9) that R=e je , i.e., \Ii\ = l. The angle 0=argument B 
is given by 

tP Q P+Q "I 

Stan-KVW+Stan-KS/TO-S taii" 1 ^/*,) 
w = l w = l n=2 _} 

where ai=j8. But since p n =-yJ(nT/b) 2 —k 2 and 8=-yJk 2 —(7r/a) 2 , one can write 

tan" 1 tfP»=wr l {(b8/ir)/(n 2 -b 2 /a 2 y /2 }. 

After similarly expressing the other arctans in terms of arc sines, one can rewrite (11) as 

-=-+23 sin 1 , —s, sin * , >; sin 1 , ' ' =. (\o\ 

2 2 n ^ 2 V^ 2 -l »=i ->Jn 2 -(b/a) 2 n=i -y/n 2 -(c/a) 2 K } 

A similar expression for the exact 6 ex has been derived from (10) and is given in the following: 
% ' «« (i_* ln « « i n «W f!»j i, qW, (™; *, oVs, (*; £, o) 

2 2 7r \ a b a cj \w / \w a J \ir a / 

SAu; v, o)= J- |_sin- (n2 _ i)2J1/2 --j (13) 



where 



We shall now go on to show that the expression for the reflection coefficient derived from 
the truncated set of equations becomes identical with the exact solution as P— ><», Q-^°°, 
provided we choose P/Q=b/c. To this end, consider (12) which gives the angle of the reflection 
coefficient 8 in terms of the arc sine series. Rewrite the expression for 6 in (12) as 

2-2 + \^ 2 sm \j?=i) Tt^nf \£i Sm VV^^W *£inf 

f Q . _, cS/tt c5 Q n aS.oS^l 65^1 c8Al ^ 

-J Ssin x ' x, — ~S- X--+- S --~S---S"' (14) 

^n=i ■yjn 2 —(c/a) 2 if n=inj ir tt n=i n it n =i n w n =in 

Now, use the well known relation [Magnus and Oberhetinger, 1954], 

lim ( ^ — InM )=7=Euler , s Constant 

in (14) after letting P-^oo, Q->cx>, and derive 

Km *=*°*+S 2 (8a/T', l,0)-S l (8b/Tr; b/a, 0)-S 1 (8c/t] c/a, 0) 
p-»<» 2t 2i T 



If we now let P/Q=b/c, we can readily derive 



+lim/- ln (P+Q)—- ln P-- ln Q V (15) 

P_>oo ^ 7T 7T 7T I 



lim e/2=~ad/ir |~1-- In ~- ln -l+S^ta/ir; 1, 0)-S,(«&/t; 6/o, 0)-S l (de/ir,c/a, 0). (16) 
p_>oo 2 |_ a o a cj 



q^«>,p/q=£ 
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Comparison shows that the limiting value of B agrees identically with (13), which is the 
exact expression given by Hurd [1954]. 

It is also interesting to note that if P and Q are increased indefinitely, such that C=P/Q 
does not equal b/c, the expression for 6 approaches the limit 

lim J=4~ [l-~ In (l +£)-*- In (l + C)~]+S 2 (5a/T] l,0)-Si(lb/r', b/a, 0)-S 1 (8c/tt; c/a, 0), 

Q->» ( 17 ) 

which is different from the correct solution. Hence, even though the solution for the reflection 
coefficient will converge to a limit when larger and larger size truncated matrices are taken, 
the convergence is relative and the solution obtained will be incorrect except for the unique 
choice of the ratio P/Q=b/c. In what follows, we go on to demonstrate that this particular 
choice of P/Q ratio is also consistent with the criterion based on the edge condition at the edge 
of the bifurcation. To this end, we first study the behavior of the higher order coefficients, 
A r , in the following section. 

5. Asymptotic Behavior of Higher Order Coefficients 

In this section, starting from (8), we shall develop an asymptotic expression for A T /A 
for large P and Q. 

Let us rewrite (8) as 

A ,_ £ (1-ttr/ftQ £ (l-«r/7») P +5"> (l+«l/«n) 
7 „ = i (1+tti/AO n = l (l+ai/Tn) n-1 (l~<*rK) 

where 

^— (4/J)™~jjfe (18) 

' simrb/a 

Since 

717T UIT , 717T P ! 

p n — > -j-j T w -> — ■ and a n -> — lor large n, 
be a 

it is convenient to introduce some additional factors in the numerator and denominator of (18) 
and rearrange it as 

r \n=i 0—<Xrblnw) „-i (l—a T c/nir) w =i (1— <*rK) / 

f £ (l+ ttl ftM « (l+ ttlC /7l7T) W« (l+«l/«n) ^ 

\n=i (1+ajPn) n=i (l+ai/7n) w =i (l+axd/mr) j 



(l+ai/P n ) n=i (l+«i/7n) n=l (1+axd/mr) f 

P Q P+Q 

M , v II (l — arb/nw) II (l—a T c/7iT) U (1+a^a/mr) 

{l—a r a/r7r) n =i w= i TO=1 



(l+ ai a/nr) p +^ (l-^a/riTr) 5 (l+a^/rwr) H (1+arf/rMr) 

n=l n = l 7i=l 



(19) 



It may be shown that the factors inside the first two (braces) tend to constants as P and 
Q are increased indefinitely for a given r. This is proven in appendix 1. 

It will therefore be sufficient to study the behavior of the ratio of the products F, where 

p Q 

II (l—a r b/n7r) II (l—a r cfnir) 

F= n -=± ^ *=! - (1-o^nr), (20) 

II (1-— a r a/mr) 



n=l 



and of a similar ratio of products with —ax replacing a r . 
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Let us rewrite (20) as 

p Q 

P«rV*ii(i- ar b/n<jr)Q a r c /* n (l-arc/mr) /pi/» )ar „ /T 

(P+#) a ' a/ * n (1 -a r a/nT) V 

71=1 

and then use the following representation [Magnus and Oberhetinger, 1954] of l/r(2+l) where 
r(#) is the Gamma function of argument x. 

i m 

-—- --==limm- z n(l + z/ft), (22) 

i- V'S"! 1 / m->oo 7Z=1 

to recast (20) into its asymptotic form for large P and Q. We derive using (22) in (21), and 
letting P/Q=C (constant), 

limF=(l-« r a/rir) rn [t!^ „. , -(l + l/C)»^(l+C)^ / '. (23) 

p-^oo 1 (1 —a r 0/7r) 1 U —a r C/7r) 

Q->oo 

Next, using [Magnus and Oberhetinger, 1954] 

r(i-«)« 



sin7T2r(2;) 
in (23) one derives 

Urn *-=(l/T).(l- ar a/nr)- 8in M) Si " ^\ T ^ T ^ M .(l+UO)^(i + C)^. (24) 
p_>oo sin a r a 1 \a r (i/7r) 

The next step is to study the limit of F as r, hence a r becomes very large. We have from 
Stirling's formula, 

r (z) ~ (2tt) ^ ^i ^~ 2 for a large 0. (25) 

Using (25) in (24) and letting a T -^nr/a we find 

/>&\ r»/« (re\ rc/a 

F-^LK r ^ r^~ jk (l + l/CT /a (l + C) rc/a for large r, P, and & (26) 

where L is a constant independent of r and K T is an oscillatory function of r such that K r -f->0 
and |2? r |<l. Note that we have used the fact that { 1 — (a r a/ir)/r}/sm a r a has the limit (— ) r+1 
(1/nr) as a r -^rir/a. 
Writing 

~r_f __L?^ 1 rb/a+rc/a 

r ~\a^a) 
we obtain from (26) 

LK, (l+l/g) fft/a (l + g) rc/g r27 x 

r 3 /2 (i+c/6)^/«(i+j/ c )^Ai ^' ; 

for P/Q= C, P-> °° , Q-> °° , and r large. 

Equation (27) gives the desired asymptotic behavior of F. It is not difficult to show that 
the ratio of products 

II (l+oLxalnir) 

= G,say, 



(l+«,a/nr) £ ( 1+ai j /rMr ) n (l+a,e/rwr) 
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has the asymptotic behavior 

lim g ^i-L 1 / ■ r ( 1 +"'^ r ( 1 + g ' c /') (l+l/C7)--.»/-(l+0)-V/- (28) 

p^co (l+aia/rx) r(l+aia/7r) 

and hence tends to a nonzero constant as r becomes large. We can sum up all this and arrive 
at the asymptotic behavior 

um ^ r _u ^ /2 (1+c/6)rt/J(1 + b/c yc /a ) lor lar S e r - W) 

Q->oo 

Letting X= b/c, we can rewrite (29) after some manipulation as 

Km A'~0 (( r )~ 3/2 {(^)(f J 7 "}') for lar S e r ' ( 3 °) 

Let C=P/Q=<rX. Substituting for C and rearranging (30) one may obtain 

(^)(?) &/a==={1 + ((T ~ 1) ^ }(7 ~^^ Say ' (31) 

withp=b/a=X/(l+X). 

Using (31) in (30) we write 

lim 4;~0((r)- 3 /y). (32) 

P-*co 

Q->co 

Without loss of generality, let us let b^>c implying X=b/c>l and hence 0.5<p<O- It is 
shown in appendix 2 that under this condition 

77>1 for o-^l 

and b}^ direct substitution we see that 

77 = 1 for cr=l. 

Using these in (32), it follows that the higher order coefficients A r have an exponential growth 
for 0-3^1, i.e., for C=P/Q^X=b/c. Only for cr=l, which implies P/Q=b/c, A' r 's exhibit an 
algebraic behavior and are 0(r~ 3/2 ) for large r. 

6. The Edge Condition and the Proper Choice of C 

The condition at the edge of septum requires that the field potential <p go to zero at the 
edge as r 1/2 and that \V<p\ go to infinity as r~ 1/2 . Now it is seen from (2a) that 



<><Pa 
dz 



7T0 _°° 

= -a 1 Ae- a i z sin — hZ) i A n sin (mrb/a)} • a n e a «*. (33) 

x = b Q> n=\ 



It is obvious that if A' n = — sin (^7r&/a).4 w /sin (7r6/a) increases exponentially for large n, 
which it does when P/Q==C9 2£ X=b/c, the expression for d(p A /dz\ x=b at the edge goes to infinity 
in a much stronger manner than z~ 1 ' 2 as z^O. It is only when C=X and ^ is 0(l/?r /2 ) for 
large n that the sum of the series is 0(2~ 1/2 ) as z^O, as may be shown by following a method 
due to Hurd. It may be pointed out that in the exact solution obtained by Hurd A' n is 0(l/n 3/2 ) 
for large n, as it of course must be. 
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The conclusion is then that the edge condition is satisfied only with the choice of C=X, 
i.e., P/Q=b/c. When P and Q are increased indefinitely, consistent with the condition C=X, 
the solution converges to the one which satisfies the required behavior at the edge and thus 
yields an unique answer satisfying the physical conditions. 

7. Conclusion 

It has been shown that the solution of a doubly infinite set of equations may exhibit relative 
convergence in that it may be dependent on the manner in which the doubly infinite set is 
truncated. An additional condition, such as the edge condition, may be used to determine 
the correct way of truncating the infinite set. Only with the proper choice of the ratio of the 
number of equations, the solution of the truncated set yields the exact solution as the number 
of equations is increased indefinitely. Otherwise, the solution converges to an incorrect 
answer which depends on the ratio of the number of equations from the doubly infinite set. 

It is hoped that this work may provide insight for handling a class of problems of a similar 
type which do not lend themselves to an exact solution. 
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8. Appendix 1. Convergence of Certain Products 

We shall prove the convergence of the products of the type 

n ()-«r"l™) (A _ 1} 

n=i (l—a r /a n ) 

as P— ><». These products appear in (19) in the text. To prove the convergence of the 
products in (A-l), it will be sufficient to show that the series S(a r ) given by 

S(a r )=jt {In (1 -a r a/mr) -In (1 -a r fa n ) } (A-2) 

72 = 1 

converges. Consider the Rth terms of the series for R large such that (c^a/JiV) < <0 • Expand 
the log functions in a series and obtain 

In (l-a r a/ Pit) -In (l — a r /a B ) =a r (a/Rir—l/a R ) + 0(l/R 2 ). (A-3) 



Using a R = ^(Rir/a) 2 —k 2 one derives 

l/a R = a/R7r+0(l/R 3 ). (A-4) 

From (A-3) and (A-4), it is obvious that the i?th term of the series in (A-2) is at least 0(1 jR 2 ) 
and hence the series is convergent. By exponentiating the series, one immediately concludes 
that the product under consideration is convergent. Similar reasoning applies for the other 
products appearing in the braces in (19). 

9. Appendix 2 
To prove the inequality 

l + ((j-l)^>c7 p , 0.5<^<1 (A-5) 

253 



for a 9^1. Note for <7=1, the above becomes an equality. 

We shall take a geometric approach to show that the above is true. The two functions 
l-\-(<r—l)p and u v are plotted as functions of a in figure 2, and it is seen that they have a 
common point at <r=l. 



Figure 2. Plot of functions appearing in inequality 
relation. 




The slope of the curve a p satisfies the following: 

Slope of a p =pa p ~ 1 

<p for (7>1 
>p for cr<l 
=p for o-=l. 



(A-6) 



Since the slope of the straight line is p and the two curves have a common point at cr=l, it 
follows from above that the line l-\-(<r—l)p is tangent to the curve a p at a=l and, further- 
more, the curve a p lies below this tangent for all c^l. 
Hence, the inequality is proven. 
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